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; The call is standard call-by-reference. MTHSGSINCOS_R7 is a special 
; routine which is the same as MTHSGSINCOS, except a faster non- 
; : are saved. call oy used with the argument in *RO and no registers 
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FUNCTIONAL DESCRIPTION: 


The GSIN, GCOS and GSINCOS routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = [1*(2epi) + Ie(pi/4) + YI, 


where I1 and I are integers, 0 =< I < 8 and 0 =< Y1 < pi/4&. Since GSIN and 
GCOS have a period of 2*pi it follows that 


GSIN(x) = GSIN(I*(pi/4) + Y1)) and 


GCOS(x) = GCOS(I*(pi/4) + Y1)). 
Using the trigonometric identities for the sum and difference of two angles, 
the following table can be generated: 
If l= then GSIN(x) = and GCOS(x) = 

0 GSIN(Y1) GCOS(Y1) 

1 axes opi ( oor? GSIN(pi/4-Y1) 

¢ GCOS(Y1) -GSIN(Y1) 

ss ieigi reer l -GCOS(pi/4-¥1) 

4 -GSIN(Y1) -GCOS(Y1) 

5 “Sens eierr “"Srineree ee 

6 -GCOS(Y1) GSIN(Y1) 

7 -GSIN(pi/4-Y1) GCOS (pi/s-Y1) 
Let Y be defined as Y = Y1 if I is even and Y = pi/4 - Y1, if 1 is odd, then 


each entry of the above table is of the for +/-GSIN(Y) or +/-GCOS(Y). Based 
on the above remarks, the GSIN, GCOS and GSINCOS routines process the input 
argument x, to obtain 1 and Y, and based on | selects a suitable polynomial 
approximation, p(Y), to evaluate the desired fuction. 
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INPUT PARAMETERS: 


00000004 LONG = 

00000004 x = 1*LONG ; x is input angle in radians 
00000008 sine = Gotane ; sine is GSIN(x) 

0000000C cosine = 3*LONG ; cosine is GCOS(x) 
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; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:30:46 VAX/VMS Macro v04-00 e 
ATHSGSINCOS © Radian arguments oreee= 1b Tiidesbe EMTANTS. Saeaarncsinco.maR;1  o (3), 
E 
FD 359 
E 0; , 
: 1 ; Return sine and cosine of argument 
i 
OOFC ; 5 -ENTRY MTHSGSINCOS, “M<R2, R3, R4, RS, R6,R7> 
E 39 MTHSFLAG_JACKET 
60 00000000 ' GF 9E E MOVAB G*MTHSSJACKET_HND, (FP) 
E ; set handler address to jacket 
E9 3; handler 
£y 
E9 38 
0 04 BC 5OF E9 39 MOVG ax (AP), RO 
OOO0036E'EF 16 O2EE 340 JSB MTHSGSINCOS_R7 
8 BC 0 O2F4 41 mova RO, a@sine(AP) 
oc BC 52 7D O2F8 rt: mMOVa R2, @cosine(AP) 
O2FC 4 RET 
O2FD 44 
O2FD 45 
O2FD 46 
O2FD 47 -SBTTL MTHSGSIN 
O2FD 48 
O2FD 49 ; 
O2FD 50 ; Return sine of argument 
psed 51; 
O2FD 26 
O2FD 5 
OOFC 8 4 338 -ENTRY MTHSGSIN, “M<R2, R3, R4, RS, RE, R7> 
0 FF 56 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 8 FF MOVAB G*MTHSSJACKET_HND, (FP) 
0306 3 set handler acdress to jacket 
Babe 3 handler 
306 
3308 57 
50 04 BC SOFD 0306 58 MOVG ax(AP), RO 
00000468'EF 16 0308 59 JSB MTHSGSIN_R7 
04 0311 60 RET 
1 61 
1 66 
1 6 
1 64 -SBTTL MTHS$GCOS 
1 65 
1 66 3 
' ? ; Return cosine of argument 
Bi 
OOF C \¢ 4 -ENTRY MTHSGCOS, “M<R2, R3, R4, RS, RO, R7> 
1? rg MTHSFLAG_JACKET 


BTNSGSINCOS 
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Ronn) Se ee 


re 


6D 00000000'GF MOVAB G*MTHSSJACKET_HND, (FP) 


; pot Bendter address to jacket 


ccococcce Ss 


4 
B 
8 : ha 
8 
5 4 BC SOFD B a MOVG ax (AP) 
x 
00 oosra: EF 16 8 2 J$B8 MInSGCOS” Q, 
04 : a RET 


GTR eeMcos 
2-00 
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4-00 
SD = Degrees EP-1984 4:04 (CMTHRTL.SRCIMTHGSINCO.MAR; 1 
-SBTTL MTHSGSINCOSD - Degrees 


FUNCTIONAL DESCRIPTION: 


The GSIND, GCOSD and GSINCOSD routines are based on octant reduction. Given 
an argument, x, it is written in the form 


x = 11#360 + 1*45 + Y1, 


where 11 and 1 are integers, 0 =< I < 8 and 0 =< Y1 < 45. Since GSIND and 
GCOSD have a period of 560 it follows that 


GSIND(x) = GSIND(1*45 + Y1) and 
GCOSD(x) = GCOSD(1#45 + Y1). 


Using the trigonometric identities for the sum and difference of two angles, 
) 


the following table can be generated: 
If Il = then GSIND(x) = and GCOSD(x) = 
0 GSIND(Y1) GCOSD(Y1) 
1 GCOSD(45-Y1) GSIND(45-Y1) 
¢ GCOSD(Y1) -GSIND(Y1) 
GSIND(45-Y1) -GCOSD (45-71) 
4 -GSIND(Y1) -GCOSD(Y1) 
5 -GCOSD(45-Y1) -GSIND(45-¥1) 
6 -GCOSD(Y1) GSIND(Y1) 
7 -GSIND(45-Y1) GCOS (45-¥1) 
Let Y be defined as Y = Y1 if I is even and Y = 45 - Y1, if I is odd, then 
each entry of the above table is of the for +/-GSIN(Y) or +/-GCOS(Y). Based 


fr) 
on the above remarks, the GSIND, GCOSD and GSINCOSD routines process the input 
argument x, to obtain I and Y, and based on | selects a suitable polynomial 
approximation, p(Y), to evaluate the desired fuction. 


00000004 LONG = 4 
00000008 sind = s,i-oue 
0000000C cosd = 5*LONG 
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MTHSGSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:30:46 VAX/VMS Macro v04-00 P 
BOOS ATHSGSINCOSD = Degrees Seem 8bs P1:S208 PAYANTS Bacto VOGTOO wars 29° 
OOFC 7 84 -ENTRY MTHSGSINCOSD “M<R2, R3, R4, RS, RE, R7> 
; 425 MTHSFLAG_ JACKET 
6D Q0000000'GF 9€ 9 MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 
3; handler 
50 4 BC SOFD ¢ $ MOVG ax(AP), RO 
00000584" F 16 4 8 JSB MTHSGSINCOSD_R7 
08 BC 8 a3 B 4 MOVaQ Rg. a@sind(APy 
Oc BC 7 if ? ? MOVa R2, a@cosd(AP) 
04 t3 4 ¢ RET 
44 4 
44 434 
44 435 
OOF C rr} 2 § -ENTRY MTHSGSIND “M<R2, R3, R4, RS, RE, R7> 
0346 438 MTHSFLAG_JACKET 
6D 00000000 ' GF SE pace MOVAB G*MTHSSJACKET_HND, (FP) 
034D 3 set handler address to jacket 
034D 3; handler 
034D 
034D 439 
50 04 BC SOFD 0340 440 MOVG @Xx(AP), RO 
OOO00SE2"EF 16 0352 441 JSB MTHSGSIND_R7 
035 re 
04 035 44 RET 
0359 3444 
359 8445 
$28 £68 
OOFC 59 44 -ENTRY MTHS$GCOSD “M<R2, R3, R4, RS, RE, R7> 
0358 448 
0358 449 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 038 MOVAB G*MTHSSJACKET_HND, (FP) 
36 3 set handler address to jacket 
B30 3; handler 
36 450 
04 BC SOFD 3 451 MOVG @X(AP), RO 
00000648'EF 16 036 43¢ JSB MTHSGCOSD_R7 
36D 386 45 
04 36D 454 RET 
36E 455 
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O-008 ATHSGSINCOS. R7 See 1982 te 4:04 MIMRTL Se RCM MTHGSINCO.MAR; 1 . (10). 
6 ish .SBTTL MTHSGSINCOS_R7 
4: 4 : ; This routine computes the GSIN and GCOS of the G-format value of RO/R1. The 
6€ 289 3 gomputet ton is performed one of three ways depending on the size of the | 
43 ie ; input argument, X: 
; 
6—€ 465; 1) If iki < pi/4, then X is used directly in polynomial approximation 
6E 464: of GSIN and GCos. dipdice ” 
6— 465 ; 2) Tt pi/4 =< ixt < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
6€ £96 3 to reduce the argument to an equivalent orgynens in radians, Ye one 
6— 467 ; the octant, I ceatetning the Srouaent. is then evaluated i 
6E $98 3 RF ene noavets chosen as a function of I, ae, compute GSIN(X) and GCOS (x). 
6E rt 3 3) =< iki, then the Sareut ial REDUCE_LARGE is called to 
8 6— 470; legesk te argument to an equivalent argument in cycles, Y ‘a the 
6— 471; octant, I, contanining the argument. is then evaluated 
Baee tis ; polynomials chosen as a function of 1, to compute GSIN(X) ne" GCOS (x). 
O3ee 474 MTHSGSINCOS_R7:: 
56 50 SOFD 395 475 MOVG RO, R6 ; R6 = X 
10° 18 03 476 BGEQ POS _SINCOS ; 
0000037F "EF 16 0374 477 JSB SINTOS 3 RO/R1 = GSINCiX!), R2/R3 = GCOS(X) 
50 50 52FD O37A 478 MNEGG RO, RO ; RO/R1 = GSIN(X) 
05 O37E 479 RSB 
37F $30 
7F «481 SINCOS: 
50 8000 8F AA 7F tee #*x8000, RO ; RO/R1 = ix! 
84 483 POS_SINCOS: 
50 FC7F CF SIFD Baar ret CMPG G_PI_OV : Compare pi/4 with iX: 
14 038A 485 BGTR shat SThéoS ; No argument reduction is necessary 
50 FC7F CF SIFD O38C 486 CMPG _PT_Ov_ 3; Compare 9*pi/4 with {Xi 
3 618 0392 487 BGEQ 
OOAE 31 0394 488 BRW 4 pint ; Use special logic for ix! > 9=pi/4 
0397 489 
839) 490 ; 
397 491 =) pi/& =< IK! < Depi/4 
Ba3e $36 : 
00000685"EF 16 0397 493 is: JSB REDUCE _MEDIUM ; Medium argument reduction routine 
0390 494 3; R4/R7 = Y = reduced argument 
0390 495 3; R2 = octant 
7E 54 7DFD O039D 496 MOVO R4, =(SP) ; Save reduced argument on stack 
52 OD Al 497 PUSHL R2 3; Save octant bits on stack 
sasoant | att 16 AS 498 JSB M_COS 3 te = GCOS(X) 
52 E00 Ao 499 MOVL (SP)+, R2 ; R2 =_Octant bits 
4 7DFD A 0 MOVO (SP)+, R4 ; R4/R7 = reduced argument 
£ 7D 1 mova = RO “pt : Save GCOS(X) on stack 
00000496" EF 16 : JSB mw §1 : RO/R1 = GSIN(X) | 
52 —E 7D B Ova (SP)+, R ; R2/R3 = GCOS(X) 
05 ; : RSB 
4 $ ; Logic for small arguments. iXi < pi/4. 
D 
BD j SMALL_SINCOS: 
50 4000 8F B61 Bb 10 CMPW #*x4000, RO ; Compare 1/2 with {Xi 
19 C 11 BLSS ; Sufficent jovecgh n9 not available 
50 3660 ? C4 6 CMPW #*X3E60, RO ; Compare with 
8 C9 1 BGEQ 1$ ; No polynomial act lentien is needed 
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6 7D 

54 g6 é 45FD 

3 4 7D 

FCDB CF 2 4 55FD 
4 f 7D 

E 7D 

FDGE CF 4 55FD 
6 44FD 

6 40FD 

BE 7 

0 

52.08 50r0 

ae. ee, 

54 0 S6FD 

4 & rt 4 

2 4 76FD 

E 2 70 

43 4 7D 

FCE1 CF 2 SSFD 
4 3 7 

4 86 

54 28 gore 

50 0 56FD 

50 54 62FD 

23 50 76FD 

2 6 7D 

-— FF 7D 

FCFF CF 98 52 S5FD 
5 BE 7D 

50 8=66E 44FD 

50 s«BE sor 
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Cosine and Sincos 16-SEP-1984 


6-SEP=1984 9 
RO, R6 
R6, Rb, RS 
R4, =(§P) 


R4 
(SP) 


RO, (SP) 
R4, #SINLENR=-1, SINTBR 
R6, R 


R6, R 
(SP)+, R2 


#1.0, R2 


R 
R2, -(SP) 


SP) 
rit: nang SINTBR 
+, R 


A R4/R5 = Xe 
she | aaah COSTBR1 


=(SP) 
#COSLENR2-1, COSTBR2 
d+, RG 
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R6/R7 = {Xi 


J 
a 
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AS I 
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R4/R7 = reduced argument (in cycles) 


R2 = octant bits 
Save octant bits on stack 


Save reduced or gument on stack 


RO/R1 = GCOS(X) 
Reduced arguaent in R3/R6 
Re aynctan atit® 


/R3 = GCOS(X) 
RO/R1 = GSIN(X) 
R2/R3 = GCOS(X) 


x*2 on stack 


50 

10 

gece et 
50 0 


50 
50 
50 


8000 8F 
FB8&3 CF 
FB8S CF 

38 


00000685 ‘EF 


07 00 52 


; Floatin 
ATHSGSIN AZ 
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; This routine computes the GSIN of the G-format value of RO/R1. 


eo 


ete Mac 
MTHRTL.SRO 
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The 


; computation is performed one of three ways depending on the size of the 
; input argument, xX 


of I to compute GS 


Y 
Y is then evaluated in a 
I to compute GSIN(X). 


1) af AY < pi/4, then X is used directly in a polynomial approximation 


0 N. 
If pi/4 =< ixi < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
to reduce the argument to an equivalent ergynent in rad 


ans, Y, and 
s then evaluated in a 


N(X). 
LARGE is called to 
in cycles, and the 


Check the sign of RO 
RO/R1 = GSIN(iX!) 
RO/R1 = GSIN(X) 
RO/R1 = {Xi 


Compare pi/4 with ‘Xi 


No argument reduction 
Compare 9*pi/4 with ix! 


Use specia 


is necessary 
logic for iXi > 9*pi/4 


Medium argument reduction routine 


R4/R7 = 
R2 = octant 


= reduced argument 


Branch to one of four polynomial 
evaluations depending on the 
octant bits. 


; 2) 
3 the octant, I, containing the argument. 
; nomial chosen as a function 
3 3) *pi/4 =< iXi, then the subroutine REDUC 
; reduce the argument to an equivalent argumen 
: octant, I, contanining the argument. 
3 polynomial chosen as a function of 
MTHS$GSIN_R7:: 

TSTG RO ; 

BGEQ POS_SIN 3 

JSB N : 

MNEGG RO, RO 3 

SB 
SIN: 

BICW #*x8000, RO 3 
POS_SIN: 

CMPG G PI_OV_4, RO : 

BGTR SMALC_SIN ; 

CMPG G_9 PT_OV_4, RO : 

BLSS LARGE _SIN i 
3 pi/4 =< Kt < 9epis4 

JSB REDUCE _MEDIUM 3 
M_SIN: CASEB R2, #0, #7 ; 
1$: «WORD P_SIN_R-1$ 3 

«WORD P_COS_R-1$ 

«WORD P_COS_R-1$ 

«WOR N_SIN_R-1$ 

«WORD N_SIN_R-1$ 

«WORD N_COS_R-1$ 

«WORD N_COS_R-1$ 

-WORD P_SIN_R-1$ 


; Logic for small arguments. 


SMALL_SIN: 
BGE 


MOVQ 


CMPW 
Q 


#*X3E60, RO 
1$ 
RO, R6 


iM! < pi/4. 


Compare with 2°-27 | ; 
No pol ynamtal evaluation is needed 
R6/R7 = X 
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O008 ATHSGSIN.R7 ety 94:30:66 MTHRTL.SRCJMTHGSINCO.MAR; 1 . (11). 
0 $0 44rd ouBy 638 MULG2 —RO,RO ; RO/RI = Xex 
FC79 CF SSFD 0488 6 POLY RO, #SINLENR-1, SINTBR ; RO/R1 = q(x*2) 
6 44FD O4BF 624 MULG R6, RO 3 RO/R1 = X*q(x*2) 
6 “or? rts 625 ADDG R6, RO 3; RO/R1 = GSIN(X) 
05 Qa? 6¢6 18: RSB 
4C 6 
4C 6 3 
4CB 629 LARGE_SIN: 
OOOOO700'EF 16 O04C¢ 650 JSB REDUCE _LARGE ; R4/R7 = reduced argument (in cycles) 
4CE 631 ; R2 = octant bit 
54 0 4C 6 ¢ L_SIN: TSTL R4 ; Check for degenerate case 
14 1 rt 4 ° i BEQL DEGENERATE_CASE_SIN 
07 00 52 8F 04D 635 CASEB R2, #0, #7 
406 6 $ 
06AA' rt | $3 1$ «WORD P_SIN_C-1$ 
0613: 408 638 -WORD P7COS"C-1$ 
0613' 040A 639 WORD P_COS"C-1$ 
O6AA' 040C 640 «WORD P_SIN_C-1$ 
06A5' et 641 «WORD N_SIN_C-1$ 
0658' 04E 648 “WORD N7COS7C-1$ 
B638, Beee 64 -WORD N_COS"C-1$ 
6A5' 0464 644 WORD N_SIN_C-1$ 
04E6 645 
04E6 646 
04E6 647 DEGENERATE_CASE_SIN: 
04E6 648 2 
52 01 8A O04E6 649 BICB #1, R2 ; Compute index as (R2 - 1)/2 
52 52 FF 8F 9C 0469 650 ROTL #-1, R2, R2 
03 00 5 8F Gers 91 CASEB R2, #0, #3 
O7AC* O4F g3¢ 1$: -WORD P_ONE-1$ 
07B8' 04F4 654 -WORD UNFL -1$ 
07B1' 0476 655 ~WORD N_ONE-1$ 
07B8' 0478 656 -WORD UNFL -1$ 
O4FA 657 
| 
| 
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4FA 659 
4FA 660 
78 96) -SBTTL MTH$GCOS_R7 
4FA 88 ; This routine computes the GCOS of the G-format value of RO/R1. The | 
4FA 664 ; computation is performed one of three ways depending on the size of the 
Te 90? ; input argument, X. The processing is the same as described for MTH$GSIN_R4. 
4FA 69 ‘ | 
4FA 668 MTHSGCOS_R7:: | 
50 53FD O4FA 669 TSTG RO ; Check for reserved operand 
50 8000 8F AA O4F 670 BICW #*x8000, RO ; RO/R1 = ix! 
50 ~—s« FBO1 gf 51FD 050 671 CMPG G_PI_ Ov_4, RO 3; Compare pi/4 with {X! 
2 14 050 ore BGTR SMALC_COs ; No argument reduction is necessary 
50 FBO1 CF S1FD OS5SOA 67 CMPG G_9 PT_Ov_4, RO ; Compare 9*pi/4 with {Xx}! 
40 19 $31 ore BLSS LARGE_Cos ; Use special logic for {Xi > 9*pi/4 
051 676 ; 
B23! 144 3 pi/4 =< IKI < 9epi/4 
00000685'EF 16 0512 679 JSB REDUCE _MEDIUM ; Medium argument reduction routine 
0518 680 ; R4&/R7 = Y = reduced argument 
0518 681 ; R2 = octant 
07 OO 52 8F 0518 ret M_COS: CASEB R2, #0, #7 ; Branch to one of four polynomial 
051C rt ; evaluations depending on the 
0514" O51C 684 1$: «WORD P_COS_R-1$ ; octant bits. 
O5A1' OS1E 685 -WORD N_SIN_R-1$ 
O5A1' 0520 686 .WORD N7SINTR-1$ 
0558' 0522 687 «WORD N_COS_R-1$ 
0558' 0524 688 [WORD N7COS7R-1$ 
05A6' 0526 689 -WORD PTSINTR-1$ 
05a6' 0528 690 -WORD PUSIN7R-1$ 
0514" 82 A 691 «WORD P_COS_R-1$ 
52C e386 
052C 693; y 
052C 694 ; Logic for small arguments. iXi < pi/4. 
052C 695 ; 
52C 696 
52C 697 SMALL_COS: 
50 4000 8F 861 52C 698 CMPW #*x4000, RO ; Compare 1/2 with {Xi . 
07 18 0531 344 GEQ 1$ : Sufficent overghang is available 
54 50 56FD 0533 00 CVTGH RO, R4 3 R4&/R7 = IX 
9508 31 0537 701 BRW NEEDS DOUBLE ; Use special logic to obtain overhang 
50 3€6 a } Ba 3A 10 1$ CMPW #*X3E60, RO ; Compare with 2*- 
C fa: oe. BGEQ 3; No potynoa’ al evaluation is needed 
50 50 44FD 0541 ie MULG2 RO,RO 3; RO/R1 = X*Xx 
FB6C CF O07 50 55FD 0545 705 POLYG RO, #COSLENR1-1, COSTBR1; RO/R1 = GCOS(X) 
oe i 
50 08 SOFD 054D 708 2% MOVG #1.0, RO ; RO/R1 = 1.0 = GCOS(X) 
05 0551 709 RSB 
55 710 
55 711 
55 ay LARGE_COS: 
00000700'EF 16 055 71 JSB REDUCE _LARGE ; R4/R7 = reduced argument (in cycles) | 
55 714 ; R2 = octant bits 
54 D5 0558 715 L_COS: TSTL WS ; Check for degenerate case | 
| 
| 
| 
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BOOS ATHSGSINCOSD_R? 6=SEP=1984 et 4:04 CMTHRTL.SRCIMTHGSINCO.MAR; 1 9° a) 2-0( 
: ; re -SBTTL MTHSGSINCOSD_R7 
584 745 ; This routine computes the GSIND and GCOSD of the G-format value of RO/R1. 
584 744 ; The computation is performed one of two ways depending on the size of the 
2 ? 2 ; input argument, X: 
584 23 : 1) If iki < 65, then X is used directly in polynomial approximation 
584 748; of GSIND and GCOSD. 
fe 749 ; 2) If 45 =< ixt, then the subroutine REDUCE_DEGREES is called to reduce 
B2 4 750; the argument to an equivalent argument in degrees, Y, and the 
584 751; octant, I, containing the argument Y is then evaluated in two 
R284 P36 ; polynonials chosen as a function of I, to compute GSIND(X) and 
584 ; COSD(X). : 
0584 754 
Rage 755 MTHSGSINCOSD_R7:: 
50 S3FD 0584 £36 TSTG 
10 18 Q587 75 BGEQ SINCOSD 
50 8000 8F AA 0589 758 BICW *x8000, RO : RO/R1 = IX! 
O0000599"EF 16 OS8E 759 JSB SINCOSD ; RO/R1 = GSIND(IX!) 
0594 760 ; R2/R3 = GCOSD(IX!) 
50 50 52FD 0594 761 MNEGG RO, RO 3 R2/R3= -GSIND(IX!) 
05 0598 poe RSB 
0599 76 
0599 764 SINCOSD: 
50 FA92 CF SIFD 0599 765 CMPG G45, RO ; Compare 45 to {X! 
4 14 Q59F 766 BGTR SMALL_SINCOSD i special processing for small arg | 
00000974"EF 16 O5A1 767 JSB REDUCE_DEGREES ; R6/R7 = reduced argument | 
O5A7 768 ; R3 = octant 
oe 7D QSA7 769 MOVQ RG. - (SP) ; Save reduced arg 
5 DD OSAA 770 PUSHL R ; Save octant bits 
oo009Ese "eT 16 OSAC 771 JSB EVAL_COSD ; RO/R1 = GCOSD(Y) 
5 BE 00 0582 Ak: MOVL (SP)F, R3 ; R3 =_octant bits | 
56 6—€ 7D 0585 77 MOVQ (SP), R6 ; R6/R7 = reduced argument 
6— 50 7D 0588 774 mMOVa RO, (SP) 3; Save GCOSD(Y) 
00000605 ‘EF 16 0588 775 JSB EVAL_SIND ; RO/R1 = GSIND(Y) 
52s BE 7D 05¢1 776 MOVQ (SP)#, R2 ; R2/R3 = GCOSD(Y) 
05 B2te 777 RSB 
5¢5 (778 
v5C5 779 ; 
05¢5 780 SMALL_SINCOSD: 
| oe 8208 781 SUBL 16, SP ; Allocate 4 longwords on stack 
6— 50 7D 0O5C8 At MOVQ . (SP) ; Save argument 
00000672'EF 15 O5¢B 78 JSB srt cOSD ; RO/R1 = GCOSD(:X!) 
08 A 50 7D 5D1 eee MOVQ Tt ; Save GCOSD(:iX:) 
5 8E 7D 0505 785 MOVa (SP)+, ; RO/R1 = argument 
00000619'EF 16 0508 786 JSB SMALL_SIND ; RO/R1 = GSIND(X) 
es Fs Bade 444 MOVQ (SP)+, : R2/R3 = GCOSD(:X!) 
05 O5E1 788 RSB 
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MTHSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:30:46 VAX/VMS Macro v04-00 Page 20 MTH 
thos ATHSGSIND OR? 6-SEP=1984 tE 4:04 CMTHRTL.SRCIMTHGSINCO.MAR; 1 . (38) 2-0 
gE 730 -SBTTL MTHSGSIND_R7 
E 138 ; This routine computes the GSIND of the G-format value of RO/R1. The 
SE 795 ; computation is performed one of two ways depending on the size of the input 
ef ay 3 argument, X 
a3 £38 ; 1) af AY < 45, then X is used directly in polynomial approximation 
3 ) ° 
i: 798 ; 2) If 45 =< txt, then the subroutine REDUCE_DEGREES is called to reduce 
5E 799 ; the argument to an equivalent argument in degrees, Y, and the 
5E 00 ; octant, 1, containing the argument. Y is then evaluated in two 
2f 4 3 polynomials chosen as a function of 1, to compute GSIND(X). 
OSE $8 MTHSGSIND_R7:: 
50 S3FD OSE 04 TSTG RO ; RO/R1 = X 
, 10 18 OSE 05 BGEQ POS_SIND ; 
OOO00SF2°EF 16 OSE? 88 JSB NEG_SIND ; RO/R1 = GSIND(IX!) 
50 50 52FD OSED 80 MNEGG RO, RO ; RO/R1 = -GSIND(iX!) 
05 OSF1 808 RSB 
OSF 809 
O5F 810 NEG_SIND: 
50 8000 8F AA OSF 11 BICW #*x8000, RO ; RO/R1 = ix! 
O5F B12 POS_SIND: 
50 FA34 CF SIFD OSF7 81 CMPG P ; Compare 45 to {X! 
A 16 QSFD 814 BGTR SMALL_SIND 3 special processing for small arg 
Q0000974"EF 16 OSFF 815 JSB REDUCE_DEGREES ; R6/R7 = reduced argument 
0605 816 ; R3 = octant 
0605 817 
0605 818 EVAL_SIND: 
07 00 53 8F Q605 819 ASEB R3, #0, #7 
0679" 0609 820 18: <WORD P_SIN_D-1$ 
0546" 0608 821 «WORD P_COS_D-1$ 
05a6' 0600 8 é WORD P7COS"D-1$ 
0679' O60F 8 “WORD P7SIN7D-1$ 
0675" 0611 824 WORD N_SIN_D-1$ 
060A' 0613 825 -WORD N7COS"D-1$ 
060A‘ 061 B58 «WORD N_COS_D-1$ 
0675' 061 : 7 -WORD N_SIN_D-1$ 
0619 8 
Rel2 ; 9 
619 0 SMALL_SIND: 
50 FA22 CF SIFD 0619 831 CMPG G_SMALLD, RO ; Compare 180/pi*2*-27 with ix: 
06 14 O61F 3 BGTR 1$ 3 No polynomial evaluation is | 
56 20 7D 0621 ; MOVa RO, R6 3 =o necessary 
0658 31 0624 4 BRW P_SIN_D ; 
50 5368 627 5 1$: TSTG 0 ; Check for zero 
eh 62A § BEQL 3$ ; Return if RO = 0 
50 FA2F CF aut? 6 § CMPG § SMALLEST_DEG, RO ; Check for possible underflow on 
93 1 6 8 BLEQ $ 3: conversion to radians 
06 31 0634 39 BRW UNFL ; Underflow will occur on conversion 
52 50 FAIS CF 45FD 0637 40 2$: MULG3 G_CONVERT, RO, R2 3; R2/R3 = pi 7320 ~ 2°=-6) tix! 
50 0060 8F A2 06 4 41 SUBW #*x60, RO 7 RO/R1 = iki #2*-6 
50 52 40FD 064 4g ADDG2 = R2, RO ; RO/R1 = GSIND(IX!) = (pi/180)ix: 
05 0647 43 3$ RSB 


50 S3FD 

50 8000 8F AA 
50 = F9DB if st 
QOO00974"EF 16 
07 OO 53 8F 
Fe} 

620° 

061C' 

0581' 

0581" 

Beak. 

620° 

0540" 


Floa 
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ti 
OS 
64 
64 
64 
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.SBTTL MTHSGCOSD_R7 


; This routine computes the GCOSD of the G-format value of RO. 
i: is performed one of two wa 
; Details are given in the d 


MTHSGCOSD_R7:: 
TSTG 


1$: » WORD 


R 
#°x8000, RO 
G ee rb 


SMALL_COSD 
REDUCE_DEGREES 


R3, #0, #7 
P_Cos_6-1$ 
P-SIN“D-1$ 
N7SIN7D-1$ 
N7COS"D-1$ 
N7COS"D-1$ 
N7SIN7D-1$ 
P-SIN-D-1$ 
P~COS"D-1$ 


AX/VMS Mac 
MTHRTL.SRC 


Compare 45 to iX! 


o V04-00 
MTHGSINCO.MAR; 1 


The computation 


depending on the size of the input argument, X: 
ussion on MTH$GCOSD_R4. 


Check for reserved operand 
RO/R1 = IX! 


R6/R7 = reduced argument 
= octant 


—_ ais 


MTH! 
2-0( 


50 F9C9 CF S5SIFD CMPG G_SMALLD, RO 3; Compare 180/pi*2*-27 with ix! 
06 14 BGTR 1$ ; Check if polyinomial evaluation is 
56 50 MOVa RO, R6 ; necessary. 
S2F 331 BRW P_COS_D ; POLY needed 
50 08 SOFD 1$: nOvG #T, RO ; RO = 1. = GCOSD(:X!) 


; 
: 
; 
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MTHSGSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:30:46 VAX/VMS Macro V04-00 Page 
O O08 ATHSGCOSD>R? aes 94:30:66 EMTHRIL. SREIRTHGSINCO.MAR: 1 ° ih) | 
sae 
¢ -SBTTL REDUCE_MEDIUM 
6 4; 
6 5 ; This routine assumes that the absolute value of the argument, X, is in RO/R1 
685 § ; and Pet pi/4 =< IX! < 9epi/4. It_returns an h=- format reduced argument, Y, 
? 2 § 3; in R4/R7, and the octant bits in R2. 
3 5 39 3; The reduced ergunens is obtained by locating the octant that X is in through 
Og : ay 3 a binary search and then subtracting off a suitable multiple of pi/2 
685 35 
be 5 93 REDUCE_MEDIUM: 
50 F996 gf S1FD 0685 94 CMPG G_5_P1_OV_4, RO ; 
ae 6 B 95 BLEQ 5$ 7 IME >= Sepi/4 
50 F986 CF SIFD eee 36 CMPG G_3_PI1_OV_4, RO 
15 15 069 9 BLEQ 3$ 3 IME >= 3epis4 
50 50 5S6FD 0695 898 CVTGH RO, RO 
54 50 FOCA CF 63FD 0699 899 SUBH3 4H PI_OV_2, RO, R4 : R4/R7 = IK! = pi/2 
04 18 Q6A0 900 BGEQ 23 ; 
52 01 400 beng 901 MOVL #1, R2 3 Octant bits = 001 
05 O6A 44 RSB 
06A6 90 ; 
52 02 DO 06A6 904 2S: MOVL #2, R2 ; Octant bits = 010 
05 06A9 905 RSB 
O6AA 906 
50 50 56FD Q6AA 907 3$ CVTGH RO, RO 
54 50 F9CS CF 63FD Bene 908 SUBH3 H_PI, RO, R4 3 R4/R7 = {Xt = pi 
04 18 0685 909 BGE 4§ ; ; 
52 03 odO 0687 910 MOVL #3, R2 ; Octant bits = 011 
05 6BA 911 RSB 
0688 36 . 
52 04 DO 0688 9135 4$ MOVL #4, R2 ; Octant bits = 100 
05 6043 914 RSB 
68F 915 
50 F964 CF SIFD O6BF 916 5$ CMPG G_7_PI_OV_4, RO ; 
15 15 06C5 917 BLEQ 7$ 3 IK! >= 7epi/4 
50 50 5S6FD O06C7? 918 CVTGH RO, RO 
54 50 F9B8 CF 63FD 06CB 919 SUBH3 H_$_PI_OV_2, RO, R4 > RG/R7 = IX! = 3epi/2 
04 18 06D2 920 BGEQ «= 68 ; ; 
52 05 00 0604 921 MOVL #5, R2 3; Octant bits = 101 
05 D608 4 § RSB 
52 06 00 ay 924 6$ MOVL #6, R2 ; Octant bits = 110 
05 60B 925 RSB 
60C 926 
50 50 a4 6DC 927 7$ CVTGH RO, RO : 
54 50 F9B3 CF 63FD 6e0 928 SUBHS H_2 PI, RO, R4 ; RG/R7 = IX! = 2epi 
04 18 Q6E7 929 BGEQ = 8S ; : 
32 oF °® 6—9 950 MOVL #7, R2 ; Octant bits = 111 
0 6EC 931 RSB 
6—ED 9 ¢ ; 
52 D& OQ6ED 935 8S CLRL R2 ; Octant bits = 000 
05 6EF 934 RSB 
6FO 935 
O6FO0 936 
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This routine is used to reduce large arguments (iX! >= 9*pi/4) modulo pi/4. 
; It returns the reduced argument, Y, in R4/R5 in units of cycles, and returns 
the octant bits, I, in R2. 


The method of reduction is as follows: 


x*(4/pi) = sane te ays. here n is an integer and 1/2 =< f < 1 

= 2*(n=53)2(2*558f) #(4/pi) 

= (2°53*f)#(2*(n=-53) #4/pi) 

= K*C, where K = 2°53*f is an integer and C = 2*(n-53)*4/pi 
Let L = K*C modulo 8, where 0 =< L < 8, and let I = the integer(L) and 
h = fract(L), then if 1 is even Y = h, otherwise Y = I-h 


CONSTANTS: 


foot * aval val aleal eal talealeal al al ot ee et et et et ee ee 


WR 0 OONOA UE Wn  O OONOUN ES wr 


9 
9 
9 
9 
5 
4 
9 
9 
3 
9 
Gg 
4 
0000003D 9 LLINT_ WEIGHT = *X3D 3; weights exponent by 61 
00000020 9 WTTERM WEIGHT = “X20 ; weignts exponent by 32 
00000400 9 W-MAX QEIGHT = *x400 ; maximum unbiased exponent 
00000386 96 W-ADJOST = *X3B6 ; Used to locate binary point in 
44 4 2 10.32 ; MTHSAL_4_OV_PI_V table 
00000000 00000000 00000000 00004021 965 ~~ ~ .OCTA *%x4021 s 2°32 


966 


REDUCE _LARGE: 


The first step is to_obtain the location of the binary point in the represen- 
tation of C = 2*(n-53)*(4/pi) in two parts - the number of longwords from 

the start and the number of bits from the most significant bit of the next 
Longword. Also Kk = 2°53*f must be obtained. 


ROTL #-4, RO, R3 Shift exponext field 4 bits right 

SUBW #w_ADJUST, RS Unbias exp and adjust for leading 
sores. R3 = Location of binary 
poin 

Divide R3 by 32 and mull by 4 to get 
R4 = # of longwords (in bytes) to 

binary point. 


3. © te 
53 0386 BF A2 


54 53 FD 8F #-3, R3, R4 
54 FFFFFFO3 BF CA a°XFEFFFFO3, R4 


MOVAL G*MTHSAL_4_OV_PI_V, Ré R address of RTL vector entr 
DL V_PI_V, R 


2g 000008 O'GF oD a4. : = 
5 00099 O'GF ¢ AD G*MTHSAL_4_0 ; Re = address of MTH$AL_4_ OV_PI table 
5 54 4 SUBL R4, R 3 R goints to 1st quadword of interest 
53 €0 8F A BICB #*XxXEO, R3 ; R3(7:0) = # of bits within longword 
50 7FFO BF AA BICW #*X7FFO, RO ; Clear exponent field 
50 4150 BF AB BISW #*x4150, RO ; RO = 2*21ef 
29 4AFD CVTGL RO ; RO = High 21 bits of K 
nn. 3 ROTL #16, R1, RI ; R1 = Low 32 bits of K | 
9¢ 1 BGEQ 1$ ; Check for high bit of R1 set 
d6 INCL RO ; Adjust RO if R1 is negative 
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NOME AN 0 OD NO UE WIN  O ODONAUE WN" OOONOM 


WOOO OOOOOOOOOOOOOOOOWOOOOOOOVOOODO 
WOOO OOO OW 00000000 0909 09 09 INI NS SINS SNS SOOO 


> The next step is to generate an approximation to C, call it C** to be used 
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; in comput tng posiet C** will consist of the first three integer bits of 
; C and the first ¢ fraction bits of C. These bits will he obtained from a 
; constant stored in the interger array MTH$AL_4_OV_PI_V. 


; NOTE: The ASHQ, ADDL, and MULL instructions in the follow sections may 

; result in an nreger overflow trap. The overflow incurred is intentional, 
; so that the IV bit must be turned off. The IV bit is not restored until 

; after all of the necessary fraction bits hav been generated. 


MOVPSL =(SP) ; Put current PSL on stack 
BICL #*C<PSLSM_IV>, (SP) ; (SP) = current IV bit 
BICPSW #PSLSM_IV ; Clear integer overflow bit 


The necessary calculation to produce the reduced argument can require up to 
nine longwords of temporary work space. This work space will be allocated 
on the stack. The work space will be accessed though the use of three 
registers: R4, R5, and SP. For the purposes of comments the temporary work 
space will be referred to as locations [10 though 18. The stack and its 
pointers will look something Like this: 


OOCOCCOCGDOOOOCOCOOO 

—— "DOOD OODOOOOVWwn 
Be Ge Ge Ge Ge Ge Ge Gee 
bod 


PANO ODA NE WIN SO OD NA NE WIN OOO NO NE WIN 3 CO OD NOA UE WN OOO NOUE WO 


Ct. oS ae 
T1: 
T2: 
13: 
T4: 
T5: 
Té: 
17: 
T8: 


<---- SP 


‘ 
‘ 
' 
' 
' 
' 
' 
' 
' 
' 
‘ 

See eee e222 eee eeeeeecee ' 

' 

: 


K---- RS 


<---- RS 


The following code allocates the storage and sets up the pointers. 


SUBL #36, SP ; Allocate 9 longwords on the stack 
ADDL #16, SP, R4 ; R4& points to T 
ADDL #32, SP, R5 ; RS points to T8 


Get C** = €(0):C(1):C(2):C(3) in 15/78. C(O) though (C(3) are unsigned 
integers generated from the pinay representation of C. The high three bits 
of C{C) are the the first three bits to the left of the binary point of C. 
The remaining pits C(O) and C(1) though C(3) are the first 125 bits to the 
right of the binary point of C. Note that the C(i)'s are odjusted to 
compensate for their signed (rather than unsigned) interpretation in the EMUL 
instruction. Note also that the representation of C has no more than 1 
consequtive ones, so that no carry is possible from the adjustment. 


Wa SS SS SSS SSS SSS SSS SS 
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MTHSGSINCOS 
2-008 


57 
67 


67 


04 AE 04 AE 
64 


64 00 
04 AG = =04 AG 
08 AG 08 A4 


FFEA 57 FFFFFFFC 


6— 00 04 


own 
ooo 
eee 


Oumnwmwu 
Waa se 


04 A4 
08 A4 
OC A4 


; loss of significance when computi ng 
; of significance before converting 4/17 t and |. 
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ADDLS #12, R4, R7 
ASHQ R3, (R2), (R7) 


; sateatl ize loop counter. R7 points 


° 
; shi t — roper quadword so that 
aaah! WK TT 


SUBL #4, R7 > R7 Sotehe to 16 
2$: SUBL #4, R2 ; R2 porore fone ht Barty in 
ASHQ RB, (R2), (R7) : Shift au Hie ete ojtha C(n) is in 
3 ~ne 
BGEQ iter : Check os high bit of Cin) set 
INCL ead ; Bit set. Adjust t Cin 
3$: ACBL af PS R7, 2$ ; Loop until C 0) though C(3) are in 
: t though T8 


Generate the low 128 bits of the product o> L. This product is 
equivalent to multiplying K times C‘'' modulo 8. The result of the 
135 it is in T4/T7 with bits 31:29 of T4 the octant bits, and the remaining 
5 bits yee fs oar | bits of the product. The last 53 fraction bits (bits 
0 of T6 f 15) are non-valid fraction bits that will be used 
ik if nore gh it B.. bits need to be generated. 


mule tply the high order bits of K (RO) times C'* and store the result in 
EMUL RO, 4(R4), #0, (SP) : TO/T1 = KHI®C(3) 
EMUL RO, B(R4). 4(SP), 4(SP) : 1O/T2 = KHI*(C(2):€(3)) 
MULL3 RO, 12(R4S, (R4) : 14 = Low ah bits a KHI*C(1) 
DDL «= (RA), BC SPS : TO/T2 = KHI*C'' modulo 8 
multiply the low order bits of K (R1) times C*' and store the result in 
EMUL R1, 4(R4), #0, (R4) : 14/TS5 = KLO#C(3) 
EMUL R1, B(R4), 4(R4), 4(R4) : T4/T6 = KLOsEC (2): €(3)] 
EMUL RI. 12(R4S, B(R4S, B(R4): T4/T7 = KLO#CC(1):C€(2):C€(3)) 
MULL 1. (RS) : T8 = KLO®C(0 
ADDL = (R55): 12¢R4) : T4/T7 = KLO®C'* modulo 8 
3; Add KHI*C'' to KLO#C'' to get K*C''. Store the result in 14/17 
ADDL (SP), 4(R4) ; 
ADWC «= & (SPS, B(R4) : 
ADWC 3 3=s-: (SP). 12(R4) t 14/17 = K*C'* modulo 8 


At this point there may or may not be enough | valid bits in R3/R4 to ecerate 
Y. If the first 12 fratction atts are all 1's or O's, there a poest ility of 


Consequently. we must check for loss 


ADDL3 #*x8000, -4(R5) 
BITL #*x3FFFOOOO, 


(R5) ; If the first 14 fraction bits are 1's 
(R5) ; and the reduced arg = 1-f or the 


oo 


| 
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2-00 Repuce LARGE “SEP-1984 11:24:04 CMTHRTL.SRCIMTHGSINCO.MAR; “ta (19) 
37 «12 «O7CA 1118 BNEQ CONVERT 3 first 14 oite are 0 and the reduced 
7cc so111 3 258 = f, then (and only then) bits 
7CC «(1114 : 9: is are 5 and significance whit 
7CC «1115 ; be lost. 
re 
7¢C 1113 : More bits need to de generated to cover the loss of significance. There are 
7CC 1119 ; not enough registers to hold all the potentisi entra bits, so that the bits 
oe 1 9 3 already generated must be put on the stack. 
7CC 1122 ° 
OOOOOBSCA'EF 16 O7CC 11 ‘ JSB GEN_MORE _BITS : Generate 85 gd i tional bits and add 
702 1124 ;  ~=them to one esio ng bits. Results are 
702 1125 i; stored i 
54 04 C2 O7D2 11 § SUBL #4, RG ; Adjust R4 to A the addition of 
705 «(11 ; _ another longword of K*C' 
15 FC AS 10 €O OQO7D5 11 8 BBS #29, -4(R5), 4 ; Check if loss of significance is due 
hy 1 3 i to leading ones or zeros 
ypa 13 1 ; Lost significance due to leading zeros 
65 10A4 OF 00 EA OQO7DA 11 g FFS #0, #15, 16(R4), (RS) 3; If at least one bit is set. This 
21 12 43 1134 BNEQ co : means lost significance was minor. 
OC AG OOEFFFFFF BF D1 gre 1135 CMPL @*XEFFFFFF, 12(R4) ; If one of the five high bits is set, 
17. 15 O7EA 11 § BLEQ CONVERT ; lost significance was minor. , 
009A 31 o7ec i f BRW LEADING_ZEROS G.}. 
Ore 125 ; Lost significance due to leading ones 6.3: 
65 10 A46 OF 00 €B O7EF 1141 4S: FFC #0, #15, 16(R4), (R5) ; If at least one bit is clear. This G_7" 
OC 12 O7F5 1148 BNEQ  CONVE ; _ meanslost significance was minor. 6-9 
OC AS =6FB000000 BF D1 sCO7F7? «2114 CMPL Fantzayoores 12(R4) ; If one of the five high bits is G_9 
1E O7FF 1144 BGEQU CONV ; clear, lost significance was minor. G_cl 
11 0801 1145 BRB CPADING ONES G_M' 
Bh its es 
080 1133 ; Convert bit string to double precion reduced argument. G3! 
03 1150 ° Hu 2, 
O80 M13) CONVERT: Ke 
o 1188 3 Isolate octant bits and convert low order bits to H-format. Wo. 
65 FC AS 03 aj 6 1193 EXTZV #29, #3, -4(R5), (R5) ; TB = octant bits H-P. 
Fc AS E000 000 & CA 11 6 BICL. #*x€0006000, -4¢R5) ; Clear. octant bits 
4 55 c C3 11. «115 SUBL3 5 oe 7 R4 RS =e sot low no bits of h 
00000910" FF i 13 1138 58 Pane TO_H ; RO/ 
0c 65 § 5 1199 BLBC (R5), : er a odd or even octant bits 
! 1196 3 Octant bits are od¢c. Reduced argument equals 1 - h. 
08 50 63FD 1 1164 SUBH3 8968 Bes Y=le-h 
52 20 AE 6 1165 MOVL = 3 sp) R2 = octant bits 
0096 «31 A 1166 BRW RESTOR ., 
D 1168 ; Octant bits are even. Reduced argument equals h | 
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O O08 REDUCE * LARGE raat 1:24:04 (CMTHRTL.SRC rar MTHGSINCO.MAR; 1 . (98) | 
sao rorp BEE IB 1s: novo a0, mint os | 
: 3 = | 
52. 20 AE 0 : 1154 MOVL 8S eSB) R2 ; R2 = octant bits 
0088 = 31 1108 BRW = «=—- RESTORE 
1174 
HE 
1199 ; At this point it has been determined that there is a major loss of 
8 1128 3; significance and the processing begins a looping phase. Each iteration of 
8 1179 ; the loop will generate additional extra bits of K*C’ until enough significant 
1180 ; bits to compute Y are available. ave 'ing this time the nine longwords 
8 1 1 ; allocated on the stack will be used as follows: 
bas 1188 3 TO/T2 Temporary storage used when generating extra bits. 
0838 1186 ; T3/T7 ~=Contains all significant bits generated so far. 
$838 1187 ; T8 Contains a counter, W, indicating the appropriate exponent 
0838 1188 ; et the last longword of fraction bits used in converting 
Bs 8 1189 ; to Y. 
0838 1190 
0838 1191 LEADING_ONES: 
0838 1136 
0838 119 ; If processing continues here it is known that the loss of significance is due 
ett 113¢ ; to a string of leading ones. 
65 3D 00 0838 1196 MOVL ML_INT_WEIGHT, (R5) ; 18 = bias for beet Longword 
bate 1138 3 of t e product K*C' 
OC AS =FFFFOOOO BF D1 oa 1199 LOOP_1: CMPL #*XFFFFOOOO, 12(R4) 3 shock for enough significant bits 
28 1A 0843 1200 BGTRU CONVERT_1 ; Enough bits, Convert to floating. 
OOOOO8CA'EF 16 0845 1201 JSB are MORE _BITS 3 te/T gentores, yee 
OC AG sOFFFFFFFF 8F- D1 pase 1 sf CMPL 12(R4) ; Check for all 1 
18 1A 0853 120 BGTRU CONVERT. 1 3; Not all 1's. Enoush precision bits 
0855 1204 ; to compute Y 
64 FC A4 7DFD 0855 1205 MOVO -4(R4), (R4) Compress representation of K*C'' 
FFD9 65 20 0400 8F 3D O85A 1206 ACBW #W_MAX_WEIGHT, #W_TERM_ WEIGHT. (R5), P_1 
0862 1207 ; Increment weighting factor. If 
0862 1208 3 weighting factor is greater than 
0862 1209 3 then no more bits need to be 
Bae 1210 3 = generated. 
65 1512 
OB 1 18 ; The wetght tne factor is greeter . than 1024. This means that the reduced 
862 1214; erquaen is either not q {as *nguieneole from 1 or too small to be represented 
862 1215 ; in Feformat (i.e. underfl Zero is returned in R4 for the reduced 
62 1 16 3 argument to signal this eeetantet Note that under these conditions the 
62 1217 ; correct function value is one of the values 0, +/-1. 
62 1 18 3; correct choice is determined by the calling program besed on the octant bits 
? ! 4) 3 returned in Ri. 
534 62 1221 CLRL ; Reduced argument is zero 
52 O08 AE O03 1D ff 64 1 ¢ EXTZV 39 +f 8(SP), R2 ; R2 = octant bits 
0056 31 086A 1 BRW RESTOR 
6D 1224 
6D 1225 
| 
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ot O08 REDUCE. LARGE aes 943068 EATHRTL. SREIMTHGSINCO.MAR: 1 ° (33) 
6D 1 CONVERT_1: 
54 04 1 $6 “ADDL #4, RS : R4 points to low bits of h 
09000910" 13 20 1 8 JSB_ CVf_LONG. T0.H : RO/RS = geuen 
54 =FE74 6 0 63FD 76 «(1 SUBH3) =RO, H g-"9 32, R4 ¢ R4/R7 = 2*W2(1 = h) 
52 1C AE O03 1D EF O87D 1 0 EXTZv #29, BS; 2B(SP), R2 ; R2 =_octant bits 
54 20 AE Ce 33 1 SUBL  32($P), R4 : =Y=1- 
3A é 1 ; BRB RESTORE 
89 1 
+24 1234 
9 1235 LEADING_ZEROS: 
aa, 
39 1238 ; If processing continues here it is known that the loss of significance is due 
889 1239 ; to a string of leading zeros. Note that it is known that the loop for 
889 1240 ; leading zeros will terminate before an underflow condition occurs so that the 
ttt } } ; loop does not include a test for underflow. 
65 3D dO 0889 1 tg MOVL WL_INT_WEIGHT, (R5) ; 18 = exp bias for Last longword 
tt ! te ; of the product K*C' 
OC AS OOOOFFFF 8F D1 O88C 1246 LOOP_O: CMPL #*XOOOOFFFF, 12(R4) ; Check enough fraction bits 
19 0894 1247 BLSS CONV ; Enough bits. Convert to floating 
OOOOOBCA'EF 16 0896 1248 JSB SEM MORE_BITS 3 T2/R7 contain K*C'' 
OC AS =—DS5 «(089C 1563 1$ TSTL 12(R4) 3; Check for all O's 
09 12 O89F 1250 BNEQ CONVERT_O ; Not all O's. Enough precision bits. 
64 FC AG 7D OBA 1251 MOVQ -4(R4), (R4) ; Compress representation of K*C'' 
65 20 AO O8A5S 1 26 ADDW #W_TERM_WEIGHT, (R5) ; Increment weighting factor. 
E 11 O8AB 125 BRB LOOP _0 
OBAA 1254 
O8AA 1255 CONVERT_O: 
54 04 CO OQOBAA 1 2$ ADDL #4, R4 3 R4 points to low bits of h 
00000910'EF _ 16 OBAD 125 JSB CVf_LONG_TO_H : RO/R3 = 2*Weh 
54 50 7DFD 0883 1258 MOVO 3 R4/R7 = 2*Weh 
52 1C AE 03 1D EF O8B7 1259 EXTZV #29, #3, 28(SP), R2 ; R2 = octant bits 
54 20 AE (C2 O08B8D 1260 SUBL  32($P), R4 : RG/R7 = Y =h 
00 11 08C1 1261 BRB RESTORE 
pate 1566 
08C \$6 
th 1264 RESTORE: 
24 AE B88 O8CS 1265 BISPSW 36(SP) 3; Restore IV bit and exit 
5E 28 CO 08C6 1696 ADOL #40, SP ; Remove mask and temporary storage 
08C9 126 3; ~— from stack 
05 08c9 1 08 RSB 
BA 
CA 1270 GEN_MORE_BITS: 
O8CA 1271 
O8CA 1 6 ; : 
path 12735 ; This subroutine generotes 85 extra fraction bits and puts them to the 
BCA 1274 ; existing bits. NOTE: This routine is always entered via a JSB instruction. 
BCA 1275 ; sonsequently, SP points to the first longword BEFORE T0, rather than T0 
BCA 1 8 ; itself. 
ach 1300 3 | 
52 04 C2 OQ8CA 1 8 SUBL #4, R2 3 Adjust ointer to get next quadword | 
CD 1280 : rom MTH$AL_4_OV_PI_V 
56 62 53 79 O8CD 1281 ASHQ = RB, (R2), RO > R72 Cin) 7 
v7” «6! D1 1282 BGEQ 1$ ; Branch if high bit is clear 
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; Logic to process unsigned values greater than 2°31 - 1 
EMUL a. R7, #0, 4(SP) 


ae 008 ie 4 0 ADDL 8($P) : TO/T1 O*C(n) 
° ; = KLO* 
08 AE O08 AE a5 0 ™ EMUL Rd ° We "B(SP), 8(SP) : . 
OC AE 50 C0 ADDL BS ° 3 TO/T2 = K*C(n) 
Od 11 BRB 
; Logic to process unsigned values less than 2°31 
1$: EMUL R7, #0, 4(SP) 3; TO/T1 KLO*#C(n) 


EMUL eR: R7, 8($P), B(SP)  : TO/T2 = KeC(n) 
; Add new bits to old 
64 08 AE 


~m 
Sad 
ie 
~J 
oS 
— 

co 
~ 
mw 
wv 
~ 


0a 
CELL 
04 AE 00 27 3} fa 
08 AE 08 AE ? 0 7A 
" 
0 
0 
0 
0 
0 
0 
0 
0 
0 
" 


D3 «1283 
D3 1284 
DS 1285 
DS 1 § 
D9 1 
DD 1 8 
E4 1 
EB 1290 
BEA 1291 
EA 1 3 
EA 129 
EA 1294 
aro 1295 
BF7 1 38 
8F7 129 
co ORF? 1399 (Ra) 
04 a4 = 0€ AE D8 O8FB 1300 ADWC =. 12 (SPS, 4(R4) ; 
0 1€ 0900 1301 BCC 3$ ; Check for carry from previous add 
AS 6 ane 1306 INCL i ; Propagate carry 
03 1€ 090 138 BCC $ ; Check for carry from previous add 
AS D6 0907 1304 INCL 12(R4) ; Propagate carry 
AE DO O90A 1305 3$ MOVL 4(SP), -4(R4) 3; Move new low order bits to end of 
90F 1306 ; of old low order bits 
05 90F 1307 RSB ; 
910 1308 
0910 1309 
0910 1319 
HEE 
0910 1318 CVT_LONG_TO_H: 
Bb a 
+4 1316 : This routine converts an orrar.e of three honquerd pointed to by R4 H format. 
91 Wei 3 ; The result is returned in RO/ NOTE s routine is always entered via 
0910 1318 : a JSB unecruct len, Consequently, SP et. to the first longword BEFORE 10, 
8318 1333 3 3 rather than T0 itself. 
04 C1 $910 1331 | ADDL3) #4, SP, an? 3 RE peinte to T 
84 6EFD Ogi 1 33 CVTLH (Rd)+, ; 3 = low 33° bits of h 
07 «#18 0918 13 BGEQ 3 Check for signed conversion error 
64 D6 O9IA 3 4 INCL (Ra) ; _.and adjust accordingly 
03 1E O91C 1325 BCC 1$ ; If necessary, 
AG bg aie 1 § INCL 4(R4) 3 prop opete 
6 a2 0921 1397 1$:  SUBW #W_TERM_WEIGHT, RO i no/Rs"= Clow 43" bits of of h)/2*32 
4 oer? 924 1 3 CVTILH = (R&)+, TRO) : = 2nd lowe 2b 
eis 928 1 GEQ 2$ ; theck for signed conversion. Sad 
64 06 O92A 1330 INCL (R4) ; an ad ust accordi oF y 
66 60FD 09 1331 2$ ADDH2 (R6), RO ; RO/R low 64 bits YOR Fi 
20 A2 0930 1 : SUBW = #W_TERM_WEIGHT, RO ; RO/RS = (low 64 bit 
4 oer) 9 1 CVTILH (R&)+, TRO) ; T0/T3 = 3rd lowest 2° bits of * 
06 1 9 1334 BGEQ 3$ : Check for signed conyers ten error 
2 CF 60FD 0939 1335 ADOH H_2_T0 je (R6) $ an adjust according! y 
66 60F 93F 1 § 3$ ADOH (R67, R ; RO/R low 96 bits of h)/2°64 
0 boca 1338 soe 
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BOOS REDUCE. DEGREES SEP-1984 11:24:06 EMTHRTL. SRCORTHGSINCO.MAR; 1 °° (30) | Tab! 
ate } 9 -SBTTL REDUCE_DEGREES 
944 1 g 3 i; This routine assumes that the abeotyte value of the argument is in ars 
944 134 ; The reduction process is performed in two stages. The first stage 
944 1344 ; the reduction reduces the argument modulo 360 to a value less that 3f at 
ate ! $2 ; and the second stage reduces the argument modulo 45 to a value less than 45. 
| 
944 1 23 ; Constants used in this reduction: 
get 18 | 
944 1350 POWER_MOD 360_0: ; Powers of 2 modulo 360 for t1 = 0 | 
008 9004 9002 0001 0944 1351 ia’ paki o eprint, Vcngnae y “ | 
0080 0040 00 8 94C 1 § «WORD 16, : ° 64; 12 
OOF8 0130 0098 010 Baee } i? -WORD 256, 152, 304, 248 
48 1355 POWER_MOD_360_1: ; Powers of 2 modulo 360 for t1 <> 0 
0008 0088 0110 0088 O95C¢ 1 28 -QORD™” 136, 272, 184, 8 
0080 0040 0020 0010 0964 135 «WORD 16, 32, 64, 128 
OOF8 0130 0098 0100 Baer 1 28 «WORD 256, 182, 304, 248 
0974 135 
O974 1361 | 
0974 1389 REDUCE_DEGREES: 
50 4360 8F B1 0974 136 CMPW #*X4360, RO ; Compare ix! with 2°53 
14 0979 1364 BGTR LAST “STEP ; Branch to special logic for med arg 
0998 1366 | 
B38 tt : It is assumed here that the argument is greater than 2°53. 
0978 1 $9 : The ar ument_is reduced as follows: 
097B 1370 ; x = 2*ttf, where t 1/2 = And let 2*53*f = 
0978 1371 ; 2830841 + J2 and . = bac aE Sinc °, 2 t30 + a $4 modulo 360, we have that 
097B 1 i ; i + 3 mod ylo 360. w let Los 3° t2. Note that 
097B 1373 ; 12913)" (2°9) yeC2 wa } BE (299) 2(353) = Hence, if tl is | 
097B 1374 ; not zero K = -*5, = gitar as | Hh (2*1 *. (ee yt : Ses24t8 modulo 360. | 
097B 1375 ; For tl=0OK ee SF Consequently, define k' ofoneruens to 2*t2 if tl = | 
0978 1376 ; age Eongruent to 156*2*t2 otherwise, where 0 =< 360. Then x’ = 
bare 1 a : 2) is congruent to s modulo 360 and ‘ < 2453, 
52. 50 00 0978 1379 MOVL RO, _R2 ; R2 = high Longword of x | 
0 Sores Br Sh Seek tae BIEL SATS 49 | Aga Dg bee or 
5250 C2 (098A 1388 SUBL RO, R2- ; R2 = t'e2*7 | 
53. 397 898 1384 mova = RO, - RB ; RB/RG = J | 
a TTT rF ob u2kd 0997 1386 eusG2 RO RS a ae 
° * = 
50 9180 af A2 44 1 $ SUBW eh st RO ; RO/R1 = 64*J)1 
0 3 40FD oy ! . ADDG2 = R3, RO ; RO/R1 = 64*/1 + J2 = J modulo 45 
7. oF G 9A4 1 35 ROTL #-4, . - Be ; R2 = t' 
SRR Be ie RS SM AE ibe 
3 8S AS ze 1398 SUBW : R 2 ts 
53 BS 0983 1395 TSTW R3 ; Check for tl = 0 and choose K' 
08 12 0985 1396 BNEQ oF 3 = accordingly 
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POWER MOD_360_1CR2], R2 
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R2/R3 = K' 


R2/R3 = K' 
RO/R1 = X* (mod 45) 0 =< RO < 2°53 


: Argument reduction scheme for arguments with absolute value less than 2°53 


d erquerat Y is computed as follows: 
int(x/45) 


45 
else Y = (]4#1)*#4 


PRP ES MO 


; Compare 2°36 with ix! 
t RO/R7 = 1X!/45 


: Turn off IV to avoid an exception in EMODG 
; Move PSL to R2 


R2 
#*C<PSLSM_IV>, R2 
SLSM_IV 


pelelejelelejlelelelelelelelelelelelelelealol aleleai aol a) 


#1, #0, R6, R3, R4 


: Save current IV bit 
; Turn off integer overflow trap 


; R3 = low 32 integer bits of iXi/ 
; R4/RS = fractional part of {xi/4 


; Restore IV bit 
; R6/R7 = Integer part of {Xi/45 = I 


AUF WN (OOO NANE WN O DOONAN EWN 0 ONO Wit OOmN 


; R6/R7 = 1 + 1 


G_1_0V_45, #X_1_0V_45, RO, R3, R4 
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R3 = I = integer part of ‘{Xi/45 
arenes 5¥ ectant its are even 


R6/R7 = 1 + 1 
R6/R7 = 45*(1+1) 
R6/R7 = Y 


Save only Last three octant bits 


R6/R7 = I 
R6/R7 = -45*] 
R6/R7 = Y 


Save only last three octant bits 


Wma 
2-00 


FOAA CF 


F645 CF 


F666 CF 


F601 gf 


8000 8F 
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: Polynomial evaluation for GCOS(Y) for 


P_COS_R: 
BICw  #*xB000, R4 
CMPW = #*X4000, R4 
BGTR LESS_THAN_HALF 
NEEDS_DOUBLE : 
MULH2 Ro, RG 


MOova 
NEEDS_DOUBLE_SINCOS: 
CVTHG = RG 


POLYG R4. #COSLENR2=1, COSTBR2 
MOV (SP)+, RG 
DECW RG 


RSB 
LESS_THAN_HALF : 

CUTHG 

MULG2 R6, R6 


POLYG R6, #COSLENR1-1, COSTBR1: RO/ 


; Polynomial evaluation for -GCOS(Y) 

N_COS_R: 
BICW #*x8000, R4 j 
CMPW #*x4000, R4 : 
BGTR 1$ : 
MULH2 R4, R4 : 
mMOVa R4, -(SP) : 
CVTHG R4, RS : 
POLYG R4, #COSLENR2-1, COSTBR2; 
MOVQ (SP)+, RG : 
DECW R4 : 
SUBH2 #1, R4 : 
CVTGH RO, RO 3 
SUBH2 RO, RG : 
CVTHG R4, RO 3 
RSB 

1$: CVTH R4, R6 
MULG R6, Rb : 
POLYG #COSLENR1-1, COSTBRI; 


R6, 
XORW #*x8000, RO 


Sete Se Se Se Se Se Se 


o VO 
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; Polynomials for arguments in radians 


in radians 


R4/R7 = {Y} 
Compare 1/2 with ‘¥! 
; Sufficent overhang is available 
R4/R7 = Y*2 
Save high half of y*2 
R4/R5 = Y*2 
RO/R1 = Q(Y*2) 
R4/R7 = Y*2 
R4/R7 = Y*2/2 
R4/R7 = se - 1 
RO/R3 = Q(Y*2) 
RO/R3 = GCOS(Y) 
RO/R1 = GCOS(Y) 
72 {YI 
R6/R7 = Y*2 
1 = GCOS(Y) 
R4/R7 = {Yi 
Compare 1/2 with ‘iY! 
Sufficent overhang is available 
R4/R7 = Y*2 
Save high half of Y*2 
4/R5 = Y*2 in G-format 
RO/R1 = Q(Y*2) 
R4/R7 = va 
R4/R7 = Y*2/2 
R4/R7 = -(1 = ¥*2/2) 
RO/R3 = Q(Y*2) 
R4/R7 = -GCOS(Y) 
RO/R1 = =GCOS(Y) 
R6/R7 = {Yi 
a = y¥*2 
RO/R1 = GCOS(Y) 
RO/R1 = -GCOS(Y) 
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oe 
—2> 
eo 
> 

wm 


: Polynomial evaluation for -GSIN(Y) 


SINR: 
XORW #*x8000, R4 


; Polynomial evaluation for GSIN(Y) 


DW teee8 
nw“ 
— 
2 


n 

0 

05 OABC 

ABD 

ABD 

ABD 

ABD 

ABD 

ABD 

54 8000 8F AC 


7E ¢ 7D OAC 35 P) 3; Save high half of Y 
7E 4 rere AC P) ; Save Y in G-format 
54 6E E 45FD ACc9 SP), RS 3 R4 = Y*2 
F663 CF 6 4 55FD AC taal SINTBR ; RO/R1 = P(Y*2) 
E 44FD AD R 3 RO/R1 = YeP(Y*2) 
4 f 7D AD9 R4 : RG/R7 = Y 
50 50 56FD AD ; RO/R3 = YeP(y¥*2) 
50 54 red AE ; RO/R3 = GSIN(Y) 
50 50 76FD AES ; RO/R1 = GSINCY) 
05 OAE8 


ead ah ah sh sh sh dh th dh Ab dh Ah hh oh oh eh hh eb eh eh 
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B08 EVCLE POLYNOMIAL ; Polynomiais for argu mitt Pi te 4:04 MTHRTL.SRCJMTHGSINCO.MAR; 1 = (3%) 
m3} ! 38 -SBTTL CYCLE_POLYNOMIALS ; Polynomials for arguments in cycles 
E9 1540 
ES }8c3 
3 : ‘§ ; Polynomial evaluation for GCOS(Y) for Y in cycles 
E9 1545 ° 
EQ 1946 P_cos.c: 
54 FSBA CF 71FD E9 154 CMPH H_2_OV_PI, R4 ; Compare 2/pi with ‘iY! 
9 18 EF 1368 BGEQ 1$ ; Sufficent overhang is available 
4 4 64FD F 154 MULH2 R4, RG 3; R4/R7 = v3 
E 7D F5 1550 MOVQ R4, -(SP) ; Save high half of Y*2 
54 4 76FD F 1551 CVTHG R4, RG 3 R4&/RS = ¥*2 
F6AD CF 7 4 SSFD F 155 POLYG R4, MCOSLENC2-1, COSTBC2; RO/R1 = Q(Y*2) 
4 8 7D 0803 155 MOVa (SPs RG > R4/R? = Y* 
4 A2 0806 1554 SUBW #2, RG + RG/R7 = Y*2/4 
4 0 gerd 0809 1555 SUBH2 #1, RS 3 R4/R7 = Y*2/4 = 1 
50 50 56FD OB0D 1326 CVTGH RO, RO ; RO/R3 = Q(Y*2) 
50 54 62FD 0811 155 SUBH2 R4, RO ; RO/R3 = GCOS(Y) 
VTH RO, R ; RO/R1 = GCOS(Y) 
50 50 76FD 0815 1558 CVTHG 0, RO 0/R1 
. ii 
56 54 76FD OBIA 1561 1$ CVTHG R4, RO ; R6/R7 = {Yt 
56 = 56 44FD ah 1366 MULG2 R6, R6 3 R6/R7 = Y*2 
F647 CF 07 56 SSFD 08 156 POLY R6, #COSLENC1-1, COSTBC1; RO/R1 = GCOS(Y) - 1 
50 08 40FD 0829 1564 ADDG #1, RO ; RO/R1 = GCOS(Y) 
05 O0B2D 1565 RB 
0B2E 1298 
Be 1 
OB2E 1392 3 Polynomial evaluation for -GCOS(Y) 
OB2E 1570 ; 
BI) cos. 
54 8000 8F AA 3 1898 BICW #*x8000, R4 3 R4/R7 = {YI 
54 F570 CF 71FD 0B 1574 CMPH H_2_OV_PI, R4 : Compare 2/pi with ‘Yi 
¢? 18 0B39 1575 BGEQ 1$ : Sufficent overhang is available 
54 4 64FD pee 1376 MULH2 R4, RG 3 R4&/R7 = y*3 
7E 54 _ 7D +3 F 157 MOVQ R4, -(SP) ; Save high half of Y*2 
54 Hy 76FD 42 1378 CVTHG R4, RS : R4/RS = Y*2 
F663 CF H Pe ss sete 130) core tT ten * santa COSTBC2; ot ght = ofr 
34 A2 3 13 1 SUBW ae, R > R4/R? = Y*2/4 
4 43 1 ¢ SUBH2 #1, R4 : R4/R7 = -(1 = ¥*2/4) 
50 gro 5 15 CVTGH RO, RO ; RO/R3 = Q(Y*2) 
54 ‘2 § FD 5B 1584 SUBH2 RO, R4 3 R4/R7 = -GCOS(Y) 
50 4 76FD SF 1585 CVTHG R4, RO ; RO/R1 = -GCOS(Y) 
05 o 1 § RSB 
26 6 76FD 64 13 3 1$ CVTH R4, R6 ; R6/R7 = {Yi 
§ 6 44FD 68 15 MULG R6, R6 3 R6/R7 = Y*2 
FSFD CF 0 56 HA $f 1590 POLY R6, #COSLENC1-1, COSTBC1; RO/R1 = GCOS(Y) - 1 
50 F487 CF 50 4 fe 1591 SUBG RO, G_M1, RO ; RO/R1 = -GCOS(Y) 
0 7A 1336 RSB 
7B 159 
7B 1594 ; 
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8 
678 
78 
78 
54 8000 BF AC 
7D 
i$ 
54 4 8 
FOSC CF 55 
44 


G6 5 
Sine, Cosine and Sincos 16-SEP-1984 01:30:46 VAX/VMS M v04-00 P 
Be 913 RCI sin B 


Polynomial evaluation for -SIN(Y) 


N_C 


"yORW #*xB000, R4 ; RG/R7 = - Y 


7E 34 FD MOVO R4, =(SP) ; Save argument in H format | 
7€ 4 76FD 4 CVTH R4, =(SP) ; Save argument in G format 
a 6 FD MULG3 (S$ (SP), RG + R4/RS = Y*2 
06 4 55FD D POLYG R4 #SINLENC1, SINTBC ; RO/R1 = P(Y*2) 

BE 44FD 08 4 MULG2 (SP)+, R + RO/R1 = YeP(y*2) 
ge BS kD BOD yey. ae tsp) RENT ye 
5 BF g2FD dese SUBH2 = (SP)+, R + R4/R7 = 3/48Y 
0 6FD OBA CVTGH RO, RO + RO/R3 = Y*#P(¥*2) 
50 34 GOFD AG ADDH2 4, RO > RO/R3 = GSIN(Y) 
50 50 76FD OBAA CVTHG RO, + RO/R1 = GSIN(Y) 

° SN 
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eo 


H § 
Sine, Cosine and Sincos a jt 4 91:30:46 ee Mac 
LS - 1:26:04 


SEP=1984 


~~ oF 
rwo 

~< 
zv 
Spe 
—_. 
—> 


pene i ~SBTTL DEGREE_POLYNOMIALS 
DBAF 1 
DBAF 1 j P_COS_D: 
56 FGAG CF StF DBAF 1624 CMPG z 90_OV_PI, R6 ; Compare 90/pi with Y 
49 1 B85 1625 BGEQ g i Dou fe precision isn't needed | 
4 6 56FD 3B 16 § CVTGH R6, R4 3; R4/R7 = Y 
4 4 64FD BBB C16 MULH2 R4, R4 3 R4/R7 = ¥*2 
E 4 7D BBF 16 3 MOVa R4, -(SP) ; Save high half of Y*2 
4 54 ZoED aCe 16 CVTHG RG, RG + R4/RS = ¥*2 
F69B CF 7 4 é FD BC 1630 POLYG R4, #COSDLN1, COSDTB1 ; RO/R1 = Q(Y*%2) 
0 0 56FD BCD 1631 CVTGH RO 3 ayia = Q(Y*2) 
4 8€ 7D OBD1 1 : MOV (SP)+, RG : R4/R7 = "53 | 
4 ) A BD4 616 SUBW 8 4 3s R4/R7 = Y¥*2/2*1 
57 00001000 8F D1 OBD7? 1634 CMPL = @#* 1000, R7 : Check for loss of significanc when | 
13_ 1A +43 1635 BGTRU 1$ 3 subtracting 1 and adjust 
43 4 7DFD BEO 16 § Ovo » ~(SP) ; Save Y¥*2/2°13 on stack 
57 FFFFOOOO 8F CA BE4 16 BICL #*XFFFFOOOO, R7 ; R4/R7 = High order bits of ¥*2/2°13 
4 62FD BEB 1638 SUBH (SP) ; (SP) = Low order bits 
+3 62FD BEF 1639 SUBH (SP)+ RO : RO/R3 = Q(Y*%2) = Low pits of ¥*2/2°13 
54 62FD BFS 1640 1$ SUBH2 #1, -R4 t RG/R7 = <(1 = ¥*2/2°13) 
28 4 8 FD BF 1641 SUBH R4, RO ; RO/R3 = GCOS(Y) 
0 0 76FD BFB 1966 CVTHG RO, RO ; RO/R1 = GCOS(Y) 
OF Or) 166k site 
56 656 aaa | COO 1645 2$ MULG2 R6, R6 3 R6/R7 = ¥*2 
08 1 C04 1646 BEQL 3; Check for Y = 0 
F61B CF 07 56 S5FD cn6 1647 POLYG R6, #COSDLN2, COSDTB2 ; RO/R1 = Q(Y*2) 
= 6o oe. 8 
50 08 SOFD a¢0 1650 3$ MOVG #1, RO ; RO/R1 = GCOS(Y) 
05 C1 1651 RSB 
0C1 1636 
OCS 1884 N_coS_p | 
56 F440 CF SIFD OC13 1655 ~ ~ CMPG G_90_OV_PI, R6 3 compere 90/pi with Y 
6 2 salh $i2 1626 Pere es ia ; Double precision isn't needed 
ree OD bees 1889 MOO. ORGS tSPD poe eee 
54 4 rere C26 1660 CVTHG R4, R ; R4/R5 = Y* 

F637 CF 07 4 55FD C2A 1661 POLYG R4, #COSDLNI, COSDTB1 ; RO/R1 = Q(Y*2) 
BF ke eb ORS 1888 move” tsb) ee RG | ne See | 
ef Op he C38 1664 SUBW ’ 3 R4/R7 = ¥*2/2°13 

57 00001000 8F OD C 1665 CMPL #*x1000, R7 ; Check for loss of significanc when 
13 1A =0C4 1998 BGTRU 1$ 3 subreactig9 1 and adjust 
43 4 7DFD OC44 166 vo » ~(SP) : Save Y¥*2/2°13 on stack 
57 FFFFOOOO 8F CA C48 1908 ICL @*XFFFFOOOO, R7 ; R4/R7 = High order bits of ¥*2/2°13 
3 4 62FD OC4F 166 SUBH (SP) ; (SP)_= Low order bits 
62FD 33 1670 SUBHO §6=- (SP) +, RO : O/R3 = Q(y*2) = Low bits of ¥*2/2°13 
4 62FD OC57 1671 18 SUBH2 §6=s #1, _ RS t RG/R? = (1 = ¥*2/2413) | 
4 § FD C58 1976 SUBH RO, R4 3 R4/R7 = -COS(Y) 
0 4 ofp CSF 167 CVTHG R4, RO ; RO/R1 = -GCOS(Y) 
OF Ocee 1808 mes 
56 56 44FD C64 1676 2$ MULG2 R6, R6 3 R6/R7 = Y*2 


. 
MTHSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:30:46 VAX/VMS Macro v04-00 e 8 MTHI 
ns BEGREE. POLYNOMIALS ett OFi3 gio MTHRTL.SRCJMTHGSINCO.MAR; 1 — (de> | 1-06 
| 
D 13 0C68 1677 BEQL «6s 3$ : Check for Y = 0 
F587 cf 97 ee 55FD teh 1 28 POLYG 6, #COSDLN2, COSDTB2 =; _-RO/RI = 6COSD(Y) 
0 8000 &F ac C71 167 XORW #°k8000, RO + RO/R1 = =GCOSD(Y) 
0 i 0 RSB 
50 F384 CF SOF err : 3$: move G_M1, RO ; RO/R1 = GCOS(Y) 
C7E 1684 | 
C7E 1685 N.SIN_D 
56 56 S2FD a4 P sin pret (Me RO ; R6/R7 = -¥ 
50 (56 56 4560 CB 1 g ~""MULG3 RE, -RO, RO ; RO/R1 = Y*2 
1 ¢ 1 BEQL RE fURN 
F618 CF 9 0 S5FD 089 16 : POLYG RO, #SINDLN, SINDTB ; RO/R1 = P(Y*2) 
6 4dr (90 169 MULG2 R6, R + RO/R1 = YeP(Y*2) 
56 9 60 BF A2 0694 169¢ SUBW Ss #*K60, RO + R6/R7 = Y/2%6 
6 40FD (99 169 ADDG2 R6, RO + RO/R1 = GSIN(Y) 
0 C90 1694 RETURN: RSB 
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i § 
ine, Cosine and Sincos at tb 


- SBTTL 


MOVG 
RSB 


MOVG 
RSB 


Underflow; if 


MOVPSL 
CALLS 


BLBC 
MOVZWiL 
CLRL 


BBC 


PUSHL 
MOVZBL 


$B 
- END 


0:46 YAX/VMS M 
4:04 (CMTHRTL.S 


DEGENERATE _SOLUTIONS 


user has FU set, signal 


R2 
#0. G*MTHSSJACKET_TST 
SF S_SAVE_PSU(FP), R2 


( 
#MTHSK_FLOUNDMAT, (SP) 


a #2, G*MTHSS$SIGNAL 


° 


acro V04-00 Page 
RCJMTHGSINCO.MAR; 1 


; Answer is 1 


; Answer is -1 


error. Always return 0.0 


R2 = user's or jacket routine's PSL 
RO = TRUE if JSB from jacket routine 
branch if user did JSB 


; get user PSL saved b 


CALL 
0 = result. LIBSSIGNAL will save in 
sar’ 1 act so any handler can 
x 


up 
has user enabled Wanting under f Low? 
yes. return PC from special routine 
trap code for hardware pleat ing 
underflow convert to MTH$_FLOUNDMAT 
(32-bit VAX-11 exception Code) 


; signal (condition, PC 


return 
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aa 
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MTHSGSINCOS ; Floating Point Sine, Cosine and Sincos ots et B49 91:30:86 pAlNas Be cro V04-00 Page 41 | 
Symbol table 6-SEP=1984 11:24:04 (CMTHRTL.S RCIM THGSINCO.MAR; 1 (26) | 


SINLENC = 8? 
SINLENR = 7 
SINTBC 9 R 
SINTBR 138 R 
SMALL_CO 5 § R 
SMALL_COSD 672 R 
SMALL 4AA R 
LLL Ha R 
SMALL_SINCOSD CS R 
SMALL N $? 619 R 
FL 0 GAA R 0 
W_ADJUST = 88 
W-MAX_WEIGHT = 4 
W-TERA_WEIGHT = 000000 
= 00000004 
X_1_0V_45 = 0000C16C 
were ctewmeaeoen en oae + 
! Psect synopsis : 
PSECT name Allocation PSECT No. Attributes 
. ; 00000000 ¢ 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS 00600000 ¢ 0.) 01 ¢ 14.) NOPIC USR CON ABS LCL NOSHR EXE RD — WRT NOVEC BYTE 
_MTHSCODE OOOOOCCE ( 3278.) 02 ¢ 2.) PIC USsR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
tee er we en ee meme remeron 
! ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization $3 00 + 80 00. 06 00:00:01 18 
Command processing 109 Sb Sy Bo 3 On OP 6 
Pass 203 0:00:05.64 00:00:14.74 
Symbol table sort 0 00: eh 00:00:00.49 
Pass 2 309 B8: 88 s ¢ B28: 15-28 
Symbol table output 16 0:00:00.1 0:00:00.67 
Psect synopsis output 00:00:00.0 0:00:00.10 
Cross-reference output BB BO sth e B86; 9990 
Assembler run totals 67 00:00:10.6 0:00:32.90 | 
| 


The york ing (88 Limit was 1500 pages. 

eytes. (75 pages) of virtual memory were used to buffer gbe intermediate cod 
There were 20 pages of symbol table space gi tgcated to hold 196 non-local and 50 local symbols. | 
1790 source Lines were read in Pass 1, producin object records in Pass 2. 
10 pages of virtual memory were used to define 9 macros. | 


$uoveseesoeseesccssoesucesse$ 


Macro Library name Macros defined 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 5 
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MTHSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:30:4 AX/VMS Macro V04-00 P 
vant Macro Run Statistics 8 g-8 aatyy 9: 9:38 YRTMRTe Red THGSINCO.MAR; 1 _ (38) 


131 GETS were required to define 5 macros. 
There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGSINCO/OBJ=OBJ$:MTHGSINCO MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) #MS | 


} 
| 
' 
| 
| 
' 
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